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Nuclear structure corrections to energy levels of light muonic atoms are derived with particu¬ 
lar attention to the nuclear mass dependence. The obtained result for the 2P-2S transition of 
1.717(20) meV serves for determination of the nuclear charge radius from the spectroscopic mea¬ 
surement in muonic deuterium. 
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In order to resolve discrepancies for the proton charge 
radius , spectroscopic measurements in light muonic 
atoms, such as /xD, p, 3 He, and /x 4 He have been performed 
[4] for the comparison of nuclear charge radii with those 
obtained from traditional atomic spectroscopy or electron 
scattering from nuclei. The nuclear charge radius can be 
determined from spectroscopic measurements, provided 
the atomic structure is well known and the influence of 
nuclear excitation on atomic levels is properly accounted 
for. The atomic structure is well understood because 
one can calculate within quantum electrodynamics the 
atomic levels with very high precision, up to the value 
of fundamental constants. Much more problematic is the 
accurate description of nuclei and their electromagnetic 
interactions with surrounding electrons and muons, be¬ 
cause of the difficulty in solving quantum chromodynam- 
ics in the low energy scale. 

The nuclear polarizability effects in muonic atoms have 
been studied for some time. In 1977 Friar in Q calculated 
the nonrelativistic nuclear electric dipole polarizability 
and Coulomb corrections for muonic helium. Eighteen 
years later, Leidemann and Rosenfelder in ;(| calculated 
the inelastic contribution for /xD in a more general ap¬ 
proach by construction of the forward two-photon scat¬ 
tering amplitude for the deuteron. More recently, we 
calculated in Q nuclear structure effects in muonic deu¬ 
terium using a perturbative formalism and have shown 
the absence of the Zemach correction. The results of 
this perturbative approach have been confirmed by Friar 
in [8j using zero-range nucleon potentials. A systematic 
dispersion relation approach was used in Q to obtain 
the complete two-photon exchange contribution, but the 
result suffered from insufficient inelastic scattering data 
from the deuteron. Recently, a perturbative approach 
has been pursued by independent derivation and numer¬ 
ical calculations for /xHe in 0 and /xD in Q 

In this work we include higher order terms in the ex¬ 
pansion in a small parameter being the nuclear excita¬ 
tion energy over the muon mass, and recalculate all other 
contributions with special emphasis on the nuclear mass 
dependence and separation of the so-called pure recoil 
corrections. Since the nuclear effects are the main source 
of theoretical uncertainties in muonic atoms, we aim to 
calculate them as accurately as possible, in order to ex¬ 
tract precise nuclear charge radii from the muonic atom 
spectroscopy. Our main limitation will come from the 


simplified model of nuclear interaction with the electro¬ 
magnetic field which assumes certain commutation rela¬ 
tions, from the neglect of possible corrections to the elec¬ 
tric dipole operator and from the uncertainty regarding 
the neutron polarizability. 

In the following we derive general formulas for the nu¬ 
clear polarizability shift using various perturbative ex¬ 
pansions. We aim to improve results obtained in Refs. 
EOS EH by correcting mass dependencies and includ¬ 
ing higher order terms. Let us first introduce the no¬ 
tation used. Positions of the muon and nucleons are 
r,r a - Corresponding relative positions with respect to 
the nuclear mass center are p , p a ■ Momenta of the muon 
and nucleons are p,p a - Relative nucleon momenta are 
Q a = Pa — Pnta/M, where the total nuclear momentum 
is P with the nuclear mass M = ^ m a . The canonical 
commutation relations 

K,pi\=iSa b S ij ( 1 ) 

for relative coordinates are the following: 

\pi,9U=i(s ab -^)S ij . ( 2 ) 


We assume that the nuclear Hamiltonian is of the form 
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where m a is a proton or a neutron mass. In what follows 
we will neglect the isospin number, so we will assume 
that each nucleon is a proton or a neutron. Under this 
assumption the electromagnetic interaction is local and 
is much easier to deal with. Later on, when matrix ele¬ 
ments are calculated for the deuteron, the correct isospin 
number is assumed. This simplified treatment is because 
the full description of nuclear electromagnetic interac¬ 
tions, including separation of the center-of-mass motion 
(T3 . have not yet been presented in the literature. 
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We start derivation from the second-order Coulomb 
interaction in the nonrelativistic approximation 


SE = ((f><t>N 
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and where Hq is the nonrelativistic Hamiltonian of the 
muon with the reduced mass 
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The distance of protons from the nuclear mass center p a is 
much smaller than that of the muon p, so the dominating 
contribution comes from the electric dipole excitations 


SE = a 2 l(j)(j)N 
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where d = X) Q= i Pa- Denoting the nuclear excitation en¬ 
ergy by E, the nonrelativistic polarizability correction is 
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The nuclear excitation energy E is much larger than a 
typical muonic atomic excitation energy, thus one may 
perform the large E expansion of the muonic matrix el¬ 
ement. The corresponding formula for this expansion is 
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From this expansion, the leading electric dipole polariz¬ 
ability contribution is [S] 


while Coulomb corrections are 
Z 4 a 6 mi 
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Sc 3 E is small and thus can be neglected for light muonic 
atoms. The dipole operator d in the above is the posi¬ 
tion of protons with respect to the nuclear mass center. 
However, one may expect some corrections to d. Indeed 
the chiral effective field theory predicts various relativis¬ 
tic corrections to the electric dipole operator. We do 
not calculate them here and therefore associate a rela¬ 
tive uncertainty of 1%, which is twice the binding energy 
per nuclear mass. 

In the evaluation of further corrections we neglect 
Coulomb corrections, and so assume the on-nrass-shell 
approximation for the muon. All corrections can there¬ 
fore be expressed in terms of the two-photon forward 
scattering amplitude. Let us consider again the related 
muonic matrix element P for the nonrelativistic two- 
Coulomb exchange 


P = 


1 


I P~ Pa | (H 0 -E 0 + E) | p- p h | 


(17) 


Using the on-mass-shell approximation and subtracting 
the point Coulomb exchange it becomes 


x (e i£ ( A "^> - i). (18) 


2 m. 


This integral can easily be performed analytically, but we 
will choose another way, which will be convenient when 
relativistic corrections are included. We will calculate di¬ 
rectly the expansion coefficients in powers of E. There 
are two characteristic integration regions: k ~ V Em and 
k ~ to, where m is the muon mass. In the first integra¬ 
tion region, where k is small, one performs an exponent 
expansion in powers of k (p a — p\). The leading quadratic 
term is the electric dipole contribution 
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and it has already been accounted for in Eq. (HUH) . The 
term with the fourth power of nucleon distances is 
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The corresponding correction to energy is 


The corresponding correction to energy is 
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These parts are due to the electric monopole, dipole, and 
the quadrupole nuclear excitations, correspondingly. 

In the second integration region, where k ~ m is large, 
one performs an expansion in powers of not exactly E, 
but of the total nuclear energy E, 


E = E + 


2 M 1 


( 22 ) 


which happens to be much more appropriate. The first 
expansion term is 


P = ^ma 2 (j) 2 (0)\f a -r , b \ 3 (23) 


and the corresponding correction to energy 
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is the modified Zemach moment. The second expansion 
term for k ~ to is 
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is proportional to k 2 , so the k integral vanishes after sub¬ 
traction of singular terms. As a result, no corrections to 
the modified Zemach moment are found. 

Consider now correction due to the finite nucleon 
size. The proton and neutron charge distribution enters 
through the convolution with the Coulomb potential in 
Eq. (fol Since their charge radii are much smaller than 
that of nuclei, one can perform an expansion of the elec¬ 
tric form factors in powers of k 2 . When k ~ y/Em, the 
electric dipole polarizability Eq. Cl is modified by 
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where 5 = ]T) a r Ea da- When k ~ to the Zemach contri¬ 
bution is corrected by 
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The case a = b is considered separately as it involves 
a momentum exchange, which is of the order of the in¬ 
verse of the proton size. When a large momentum is 
exchanged, the nucleon binding energy can be neglected 
and the muon sees free nucleons. The individual Zemach 
radii and nucleon polarizabilities are combined together 
into effective Dirac-delta type interactions and are ac¬ 
counted for in 5 NjP E in Eqs. (H51) . and (HT1) . 

Consider now corrections from the two-Coulomb ex¬ 
change using the relativistic (Dirac) Hamiltonian for the 
muon. Equation (TTSl) is replaced by 
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where E k = Vk 2 + to 2 . When k ~ y/2 to E one employs 
a small k expansion. The leading term coincides with 
Eq. (fT8l) . The next term is 

2 c ik-(f a -fj) Ek 2 

(2 to- E + k 2 ) 2 
(31) 

Only the quadratic term in nuclear distances contributes, 
and after subtraction of large k singularities the corre¬ 
sponding correction to energy 
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is in agreement with the former result of Ref. [ij in the 
infinite nuclear mass limit. 

When k ~ to one can perform the Taylor expansion 
of the integrand of Eq. (1301) in powers of E. The term 
without E, 
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is exactly the same as in the nonrelativistic limit, and 
has already been accounted for. The linear in E term in 
Eq. (1M1) is 
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The corresponding nuclear matrix element can be trans¬ 
formed using Eq. (®. so correction to energy becomes 



/ 4 7r\ 2 to (4 to 2 + 3k 2 ) 

\ k 2 J 2 tojv Ek k 2 


(35) 

a recoil correction for each individual nucleon. So, for 
k ~ m the muon sees individual nucleons and these cor¬ 
rections become the sum of nucleon recoil corrections 


S' C E = 
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Using the commutation relations of Eq. © the nuclear 
matrix element is 
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and this correction can be rewritten in the form 
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The identity (1391) is an approximate one. The electric 
dipole operator does not commute with the nuclear po¬ 
tential. Neglected terms can be interpreted as if due to 
A 2 vertex with pions. 

When k ~ to the complete two-photon exchange is a 
recoil correction from individual protons [see Eq. (32) of 

Ref. 0 ], 


with subtracted muon-nucleus recoil correction to avoid 
the double counting with the so-called pure recoil correc¬ 
tion. This is because recoil corrections are by definition 
included in the Lamb shift as a QED correction for a 
point nucleus. 

The Coulomb exchange is not a complete correction; 
there are single and double transverse photon exchange 
corrections, and their calculation is more complicated. 
The main reason for this is the overlap of the nuclear 
recoil and nuclear polarizability corrections. Let us re¬ 
peat now the calculation by replacing the two-Coulomb 
exchange amplitude, Eq. by a complete two-photon 
exchange. 

When k ~ E or k ~ \/Em one can use a dipole ap¬ 
proximation, where the coupling of the nucleus to the 
electromagnetic field is — d ■ E(R), as in Eq. ©. Correc¬ 
tion to energy due to two-photon exchange in the dipole 
approximation is [TdJ 


SE 
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The leading nonrelativistic term agrees with that in 
Eq. m, while the leading relativistic correction agrees 
with Eq. ©. The higher order correction (in powers of 
E/m) is 
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The contribution beyond the previously considered 
Coulomb part, Eq. (IMll . is 


5' h E = — a 2 0(0) In — 
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where we again subtract the corresponding nuclear recoil 
correction. The In e dependence cancels out with that 
in Eq. SO]), as it should, and the sum of higher order 
corrections is 
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These are all nuclear structure corrections up to the or¬ 
der a 5 m 2 /M. In some cases, such as for the deuteron 
nucleus, where the magnetic moment is relatively large, 
higher order effects due to the second-order magnetic in¬ 
teraction a 5 m 3 /M 2 may play a role. The corresponding 
correction for the deuteron was obtained in Ref. jj], 

SuE = i^(o)(ti0 

3 V 4to p / 

x f dE\ {(/>N\sp - s n \E) 2 , (45) 

J Et V 2 m r 

but the numerical value, as was pointed out in [0, was 
in error, so we correct it here and present our updated 
value in Table I. 
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There are, in addition, contributions due to intrinsic 
elastic and inelastic two-photon exchanges with individ¬ 
ual nucleons, which include the third Zeinach moment 
and the nucleon polarizability. While the contribution 
from the pro ton is well known from studies on muonic 
hydrogen jl5j ], AE(2S) = —36.9(2.4) /reV, less is known 
about the contribution from the neutron. Following [9], 
we assume that this contribution is as large as the in¬ 
elastic part for the proton 13.5/xeV, and associate 50% 
uncertainty. Therefore, the contribution from intrinsic 
nucleon polarizabilities and elastic two-photon exchanges 
is 


S pE 


SnE 


8Z 4 „ m: 
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TO 


^- 36.9 meV, 

rH 

3 


- 8 ^ §w l^N. 13.5 meV. 


m, 
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(46) 

(47) 


The final expression for the nuclear polarizability com¬ 
bined with the elastic contribution but with subtracted 
nuclear recoil of order a 5 m 2 /M is 


A E = SoE + ScE + SrE + SqE + SfsE + SfzE 

+SmE + SpE + 5 nE + ShqE + 5zE. (48) 


and the elastic contribution for the neutron using Galster 
parametrization [16] is found to be negligible. 

Numerical results for muonic deuterium are obtained 
by using the AV18 potential [17] with the help of a dis¬ 
crete variable representation [181 ] method for solving the 
Schrodinger equation, and are presented in Table I. They 


TABLE I: Nuclear structure corrections in muonic deuterium 
for 2P-2S transition. Fundamental physical constants are 
from Ref. jjj, and r 2 = 0.8409 2 fm 2 , r 2 = —0.1161 fm 2 . 5^ 
from El includes only the logarithmic part of SciE, which 
we find here to be a good approximation. 


Correction 

Value in meV 

Eq. 

Ref. [11] 

[UJ-AV18 

SqE 

1.910 

Cl 

iST 

1.907 

SciE 

-0.255 
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SqqE 
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-0.042 

8q\E 

0.139 
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j( 2 ) 

°D1D 3 

0.139 

8Q2E 

-0.061 
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s (2) 
°Q 

-0.061 

SfsE 

0.020 
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S (2) 

°NS 

0.015 

SfzE 

-0.018 
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5 # 

-0.017 

SrE 

-0.026 
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^ 4 0 ) + 4 0) 

-0.017 

SroE 

0.004 

(1441) 
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SmE 

-0.008 

© 

$(°) 

°M 

-0.008 

5 P E 

0.043(3) 

© 


0.0135 

SnE 

0.016(8) 

© 


0.0135 

A E 

1.717(20) 



1.681(20) 


are generally in agreement with our previous calculations 


[j] with few exceptions. Differences are due to improved 
mass dependence of corrections beyond the nonrelativis- 
tic dipole term. We also corrected the magnetic contribu¬ 
tion, which previously was in error, and included higher 
order correction SroE, and most importantly the polar¬ 
izability of the neutron. We have also included, following 
Refs. ESE3 finite size corrections 5 fs and Sfz, although 
our results are slightly different. 

In comparison to Refs EM , we agree with their nu¬ 
merics, agree with the use of reduced mass in SqE, but 
disagree with their mass dependence of all other higher 
order corrections. Moreover, our formula for the finite 
size correction <5ps slightly disagree with the correspond- 
ing Sjyg due to the opposite sign for the neutron radius 
contribution. Our SrE differs from the corresponding 
<5^°y, apart from mass dependence, due to the fact that for 
the large momentum exchange k ~ m, the dipole approx¬ 
imation does not hold and we account for this in SfioE. 
Finally, we consider the elastic contribution of the proton 
structure correction to be a part of the overall nuclear 
structure correction, in contrast to Refs. |,E2- Our 
argument is the following. The nuclear structure contri¬ 
bution to high extent is given by the forward two-photon 
scattering of the nucleus. When momentum exchange is 
much larger than the nuclear binding energy, muon sees 
individual nucleons and the total scattering amplitude 
is a coherent sum of total scattering amplitudes from 
each nucleon. By total we mean the elastic and inelas¬ 
tic contributions, since this division is a pure convention. 
Therefore, both contributions should be included in the 
calculation of the Lamb shift, and we include them, for 
convenience, in the part called the nuclear structure cor¬ 
rection. 

Considering the uncertainty related to numerical eval¬ 
uation of matrix elements, Ref. E3 h as performed calcu¬ 
lations with AV18 potentials and with various orders of 
chiral effective field theory, finding 0.6% dependence on 
the potential used. Our numerical values, when the same 
formulas are used, are in perfect agreement with those of 
Ref. EH- Since we neglect possible corrections to the 
electric dipole moment, which in fact depends on the 
model potential, we do not associate additional uncer¬ 
tainty beyond that assumed for the electric dipole polar¬ 
izability. Regarding uncalculated higher order terms, the 
most significant seems to be the Coulomb correction to 
SqE in Eq. (l2lT) . which we estimate by about 0.005 meV. 
Therefore, the final uncertainty is determined by 50% of 
the estimated neutron polarizability, 1% of A E due to 
neglect of corrections to the electric dipole moment, and 
0.005 meV due to neglected higher order terms. 

Our final result for the nuclear structure correction A E 
is not in perfect agreement with that of Ref. E2 , as 
explained above, mostly due to inclusion of the proton 
elastic contribution. In spite of other small discrepancies 
with EH, the presented perturbative approach seems to 
be more efficient than the dispersion relation approach 
of Ref. [§]. If further improvements are required, the 
best way is probably by joining within the dispersion 
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relation approach, the inelastic scattering data at high No. 2012/04/A/ST2/00105. 
energies with nuclear model calculations at low energies, 
to account properly for the high energy structure of the 
deuteron. 
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